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Abstract 

The determination of the full SUSY QCD corrections to the production of squarks and 
gluinos at hadron colliders is reviewed. The inclusion of the NLO corrections stabilizes 
the theoretical predictions of the various production cross sections significantly and leads 
to sizeable enhancements of the most relevant cross sections. We discuss the phenomeno- 
logical consequences of the results on present and future experimental analyses. Finally 
we investigate the impact of the corrections on the transverse momentum and rapidity 
distributions. 
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1 Introduction 



The search for Higgs bosons and supersymmetric particles is among the most important en- 
deavors of present and future high energy physics. The novel colored particles, squarks and 
gluinos, can be searched for at the Tevatron, a.pp collider with a cm. energy of 1.8 TeV, and the 
future LHC, a pp collider with a cm. energy of 14 TeV. Until now the search at the Tevatron 
has set the most stringent bounds on their masses. At the 95% CL, gluinos have to be heavier 
than about 175 GeV, while squarks with masses below about 175 GeV have been excluded for 
eluino masses below ~ 300 GeV0>. In the i?-parity-conserving MSSM, supersymmetric particles 
can only be produced in pairs. All supersymmetric particles will decay to the lightest super- 
symmetric particle (LSP), which is most likely to be a neutralino, stable thanks to conserved 
.R-parity. Thus the final signatures for the production of supersymmetric particles will mainly 
be jets and missing transverse energy, which is carried away by the invisible neutral LSP. 
Squarks and gluinos can be produced via the processes 

pp/pp -> qq, qq, gg,qg + X . (1) 

In the following we consider the production of all squarks and gluinos except stops. We assume 
the squarks to be mass degenerate, which is a reasonable approximation for all squark flavors 
except stops, while the light quarks (u, d, s, c, b) will be treated as massless particles. The 
production of stop pairs requires the inclusion of mass splitting and mixing effects and will be 
investigated elsewhere B. In eq. (|IJ) a summation over all possible squark flavors and charge 
conjugate final states should be implicitly understood. The calculation of the LO cross sections 
has been performed a long time agoB. Since the [unphysical] scale dependence of the LO 
quantities amounts to about 50%, the determination of the NLO corrections is necessary in order 
to gain a reliable theoretical prediction, which can be used in present and future experimental 
analyses. 



2 SUSY QCD corrections 

The evaluation of the full SUSY QCD corrections splits into two pieces, the virtual corrections, 
generated by virtual particle exchanges, and the real corrections, which originate from gluon 
radiation and the corresponding crossed processes with three-particle final statesB. 



2.1 Virtual corrections 
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Figure 1: Typical diagrams of the virtual corrections. 



The one-loop virtual corrections are built up by gluon, gluino, quark and squark exchange 
contributions [see Fig. p]]. They have to be contracted with the LO matrix elements. The 
calculation of the one-loop contributions has been performed in dimensional regularization, 



leading to the extraction of ultraviolet, infrared and collinear singularities as poles in e = 
(4 — n)/2. For the chiral 75 coupling we have used the naive scheme, which is well justified 
in the present analysis at the one-loop level. We have explicitly checked that after summing 
all virtual corrections no quadratic divergences are left over, in accordance with the general 
property of supersymmetric theories. The renormalization has been performed by identifying 
the squark and gluino masses with their pole masses, and defining the strong coupling in the 
MS scheme including five light flavors in the corresponding (3 function. The massive particles, 
i.e. squarks, gluinos and top quarks, have been decoupled by subtracting their contribution at 
vanishing momentum transfer In dimensional regularization, there is a mismatch between 
the gluonic degrees of freedom [d.o.f. = n — 2] and those of the gluino [d.o.f. = 2], so that SUSY 
is explicitly broken. In order to restore SUSY in the physical observables in the massless limit, 
an additional finite counter-term is required for the renormalization of the novel qgq vertex!^. 




The real corrections originate from the radiation of a gluon in all possible ways and from 
the crossed processes by interchanging the gluon of the final state against a light quark in 
the initial state. The phase-space integration of the final-state particles has been performed 
in n = 4 — 2e dimensions, leading to the extraction of infrared and collinear singularities as 
poles in e. In order to isolate the singularities we have introduced a cutoff A in the invariant 
mass of, say, the qg pair, which separates soft and hard gluon radiation. After evaluating all 
angular integrals and adding the virtual and real corrections, the infrared singularities cancel. 
The left-over collinear singularities are universal and are absorbed in the renormalization of 
the parton densities at NLO. We defined the parton densities in the conventional MS scheme 
including five light flavors, i.e. the squark, gluino and top quark contributions are not included 
in the mass factorization. Finally we end up with an ultraviolet, infrared and collinear finite 
partonic cross section, which is independent of the cutoff for A — > 0. 

However, there is an additional class of physical singularities, which have to be regularized^. 
In the second diagram of Fig. || an intermediate qg* state is produced, before the [off-shell] 
gluino splits into a qq pair. If the gluino mass is larger than the common squark mass, and the 
partonic cm. energy is larger than the sum of the squark and gluino masses, the intermediate 
gluino can be produced on its mass-shell. Thus the real corrections to qq production contain 
a contribution of qg production. The residue of this part corresponds to qg production with 
the subsequent gluino decay g — > qq, which is already contained in the LO cross section of qg 
pair production, including all final-state cascade decays. Thus this term has to be subtracted 
in order to derive a well-defined production cross section. Analogous subtractions emerge in all 
reactions: if the gluino mass is larger than the squark mass, the contributions from g — > qq, qq 
have to be subtracted, and in the reverse case the contributions of squark decays into gluinos 
have to subtracted. 



3 Results 



The hadronic cross sections can be obtained from the partonic ones by convolution with the 
corresponding parton densities. We have performed the numerical analysis for the Tevatron and 
the LHC For the natural renormalization/factorization scale choice Q = m, where m denotes 
the average mass of the final-state SUSY particles, the SUSY QCD corrections are large and 
positive, increasing the total cross sections by 10-90% B. This is shown in Fig. |3|, where the K 
factors, defined as the ratios of the NLO and LO cross sections, are presented as a function of 
the corresponding SUSY particle mass for the Tevatron. 

We have investigated the residual scale dependence in LO and NLO, which is presented in 
Fig. f|. The inclusion of the NLO corrections reduces the LO scale dependence by a factor 3-4 
and reaches a typical level of ~ 15%, which serves as an estimate of the remaining theoretical 
uncertainty Moreover, the dependence on different sets of parton densities is rather weak 
and leads to an additional uncertainty of ~ 10% In order to quantify the effect of the NLO 
corrections on the search for squarks and gluinos at hadron colliders, we have extracted the 
SUSY particle masses corresponding to several fixed values of the production cross sections. 
These masses are increased by 10-30 GeV at the Tevatron and by 10-50 GeV at the LHC, thus 
enhancing the present bounds on the squark and gluino masses significantly @* . 

Finally we have evaluated the QCD-corrected transverse-momentum and rapidity distribu- 
tions for all different processes. As can be inferred from Fig. |5|, the modification of the nor- 
malized distributions in NLO compared to LO is less than 10% for the transverse-momentum 
distributions and negligible for the rapidity distributions. Thus it is a sufficient approximation 
to rescale the LO distributions uniformly by the K factors of the total cross sections 
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Figure 3: K factors of the different squark and gluino production cross sections at the Tevatron. 
Parton density: GRV(94) with Q = m. Top mass: m t = 175 GeV. 
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Figure 4: Scale and parton density dependence of the total gg and gq production cross sections 
at the Tevatron in LO and NLO. Parton densities: GRV(94) (solid), CTEQ3 (dashed) and 
MRS(A') (dotted); mass parameters: nig = 280 GeV, m g = 200 GeV and m t = 175 GeV. 
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Figure 5: Normalized transverse-momentum and rapidity distributions of pp — * qq + X at the 
Tevatron in LO (dotted) and NLO (solid). Parton densities: GRV('94) with Q = m; mass 
parameters: nig = 280 GeV, m g = 200 GeV and m t = 175 GeV. 
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